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We study nonperturbative interaction corrections to the thermodynamic quantities of multichan- 
nel disordered wires in the presence of the Coulomb interactions. Within the replica nonlinear 
cr-model (NLcrM) formalism, they arise from nonperturbative soliton saddle points of the NLcrM 
action. The problem is reduced to evaluating the partition function of a replicated classical one di- 
mensional Coulomb gas. The state of the latter depends on two parameters: the number of transverse 
| channels in the wire, N c h, and the dimensionless conductance, G(Lt), of a wire segment of length 

^vq , equal to the thermal diffusion length, Lt- At relatively high temperatures, G(Lt) ^ \nN c h, the 

gas is dimerized, i.e. consists of bound neutral pairs. At lower temperatures, mTV,^ > G(Lt) ^ 1, 
the pairs overlap and form a Coulomb plasma. The crossover between the two regimes occurs at 
a parametrically large conductance G(Lt) ~ lniVch, and may be studied independently from the 
perturbative effects. Specializing to the high temperature regime, we obtain the leading nonper- 
turbative correction to the wire heat capacity. Its ratio to the heat capacity for noninteracting 
electrons, Co, is SC/Co ~ N c hG 2 (Lr)e~ 
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^ 1 I. INTRODUCTION 

The interplay between disorder and electron-electron interactions in conductors influences their low-temperature 
properties in an essential wayi£. Depending on the disorder strength, the temperature and other system parameters a 
conductor may be either in the metallic or in the insulating regime. The manifestations of electron-electron interactions 
in the two regimes are quite different. In the insulating regime the charge in a given localized site is quantized in the 
units of the electron charge, and charge discreteness effects dominate the system properties 2 . In the metallic regime 
the charge in a given volume of the conductor can change continuously and charge discreteness effects are small. The 
two regimes can be distinguished by the value of the appropriately defined dimensionless conductance G, which is 
greater than unity in the metallic regime and smaller than unity in the insulating one. If the system crosses over 
from the metallic to the insulating regime due to a change in temperature or disorder strength the Coulomb blockade 
| ' effects are expected to gradually grow and become important at G ~ 1. 

Theoretically the transition between the metallic and the insulating regimes is typically approached from the 
metallic side, G ^> 1, where electron transport can be described semiclassically. Therefore the study of incipient 
Q\ , charge discreteness effects in the metallic regime is an important problem in the theory of disordered conductors. 
lO ' This problem has recently attracted much attention J^s^ ^io, 11,12,13,14,15,16,17 ^ j n ^he metallic regime 1/G may 
be used as a small expansion parameter. For G> 1 the charge discreteness effects are exponentially small in G, and 
QO ■ their analysis requires nonperturbative methods. To date quantitative studies of nonperturbative interaction effects in 
the metallic regime have been limited to granulated systems, or to systems in which the electron-electron interaction 
is spatially separated from the disorder. The present paper is devoted the study of nonperturbative effects in the 
thermodynamic properties of homogeneously disordered wires, in which electron-electron interactions and disorder 
, spatially coexist. 

The most promising technique to study this problem is the nonlinear a- model (NLcrM), either in the replica— or 
Keldys h 19 ' 20 formulation. We use Finkelstein'si^ replica formulation of the NLcrM. We show that nonperturbative 
corrections to the thermodynamic quantities of the wire depend on two parameters: the number of channels, N c h, 
in the wire, and the dimensionless conductance, G{Lt), of the wire segment of length equal to the thermal diffusion 
length, Lt- In contrast, the perturbative corrections^ are controlled by a single parameter, G{Lt)- For example, 
the leading perturbation theory correction to the heat capacity is SCpt/Go ~ 1/G(Lt), where Co is the wire heat 
capacity in the noninteracting electron approximation. 

Within the NLcrM formalism, the nonperturbative effects are described by soliton saddle points of the NLcrM action. 
The spatial extent of the solitons is given by the thermal diffusion length Lt, and their action is equal to G(Lt)- The 
nonperturbative contribution to the thermodynamic quantities is described by the partition function for a gas of these 
solitons. We map the problem onto a one dimensional replicated Coulomb gas. At high temperatures, G(Lt) ^ In N c h, 
the Coulomb gas is dimerized, i.e. consists of widely separated neutral pairs (dimers). In the temperature range 
lniV c ^ > G{Lt) ^ 1, the dimers are ionized and form a Coulomb plasma. Since the crossover between the two 
regimes occurs at a parametrically large conductance, G{Lt) ~ \nN c h, it can be studied independently from the 
perturbative effects. In this paper we specialize to the high temperature regime, leaving consideration of the crossover 
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to the low temperature one for future work. 

The paper is organized as follows. In Sec. [II] we describe the NLcrM for multichannel wires. In Sec. IIIII we obtain 
the analytic solution for the saddle points of the NLcrM action in the limit of the infinite number of channels N c h, 
and evaluate the functional integral over the fluctuations about the saddle points. In Sec. IIVI we obtain the leading 
nonperturbative correction to the thermodynamic quantities of the wire for N c h ^> 1. In Sec. |V]we summarize our 
results. 



II. NONLINEAR u-MODEL 



We consider an infinitely long disordered wire with many transverse channels, N c h 3> 1. The disorder is assumed 
to be weak, so that the elastic mean free path I satisfies the condition kpl 3> 1, where kp is the Fermi wave number. 
We consider the temperature T to be smaller than the Thouless energy for the transverse motion, Et = D/d 2 , where 
d is the transverse wire dimension, and D is the diffusion constant. In this regime the wire is described by the 
one-dimensional NLaM. 

Thermodynamic properties of the system can be extracted from the averaged over disorder realizations replicated 

partition function, (Z p ) = (Tre~ p ~r), withp being the number of replicas. We will be interested in the thermodynamic 
potential, which can be obtained using the replica trick: 

(0) = -TQnZ) = -T lim ~ 1 . (1) 

p^O p 

In the diffusive regime the replicated partition function, {Z p ), has a functional integral representation in terms of 
NLcrM, describing the low-energy physics of the problem. The derivation of the NLcrM action has become a standard 
procedu re 18 ' 21 . Therefore, below we only present its final form, suitable for the problem under consideration. The 
NLcrM action is a functional of two fields: the Q-matrix, parameterizing the diffusive degrees of freedom of electron 
motion, and electric potential V. The former is a Hcrmitian matrix in the space of replicas and Matsubara frequencies, 
whose entries are 4x4 matrices in the space T, given by the product of spin, S, and time-reversal, T, space a 21 ' 22 . 
The slowly varying in space electric potential V a is introduced to treat the the long range part of the Coulomb 
interaction in the replica a. This part of the Coulomb interaction is of particular importance for the consideration 
below. It cannot be described by the Fermi-liquid interaction constants. Since the Fermi-liquid effects in disordered 
metals have been studied by Finkelstein~& and are not essential for the phenomena discussed in this paper, we ignore 
them in order to keep the presentation more transparent. Then the NLcrM action can be written as 

(Z p ) = J V[Q,V}e- s ®- Sc , (2a) 



S Q = A^- I ,/,'TV 



^{yQf-(e + V)Q 



+Au f dTdxY,V e ?(x,T), (2b) 
S c = ~ fdTdxdx'J2 v a(x,T)K(x-x')V a {x',T), (2c) 

■> a 

where Tr denotes the trace over the replica, Matsubara and S £g> T spaces, v is the density of states per spin at the 
Fermi level, and A in the wire cross section area. The matrices i and V have the following structure in the replica 
and S <S>T spaces: e = iS ab T3d T , V = 5 ab ToV a , with t,'s defined as n — U <g) cro, where cr,, U are the Pauli matrices 
in the S and T spaces. The term Sq defined in Eq. (|2b() , represents the part of the action that describes electrons 
moving in the presence of the auxiliary fields V a , whereas Sc, defined in Eq. (|2c|) . is the bare Coulomb action. The 
kernel K(x — x') describes the inverse effective Coulomb interaction in the wire. In particular, for a homogeneous 
wire in the absence of a nearby gate its Fourier transform is K(q) = l/e 2 \n-^p. We also assume that the external 
magnetic field is absent. The action © constitutes the NLcrM. 

The Q-matrix satisfies the nonlinear constraint Q 2 = 1. It also satisfies the charge conjugation condition 21 , 



= CQ T C T , C = 5 ab 5 e 



/O -I 
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= 8 ab 5 E£ , ® (h ® H(7 2 )), 



(3) 
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where a, b and e, e' denote replica and Matsubara indices respectively, and the superscript T denotes the transposition. 
In what follows we restrict ourselves to the case of strong spin-orbit scattering. In this case the Q-matrix belongs to 
the symplectic ensemble^, and its matrix elements are unit matrices in the spin space. 

To resolve the nonlinear constraint Q 2 = II we will use the exponential parameterization of the Q-matrix, 

Q = e lW / 2 Ae- lW/2 , kf £ , = S ab S E£ ,T sgne, {W, A} = 0, W = W + , (4) 

where {A, B} denotes the anticommutator of A and B. The invariance of the Q-matrix with respect to the operation 
of charge conjugation, Eq. ([3]), and its hermiticity impose the following matrix structure on the rotation generators, 
W, in the T-space: 

Wg = ( d e ^ ) , d+ = d, c T = -c. (5) 

The fields (d, c) ££ , represent the diffuson and cooperon degrees of freedom respectively, each being a unit matrix in 
the spin space. 

The action in Eq. (|2J) is characterized by two parameters. The first is G(Lt) = 47r TivDA/Lt, where Lt = \JTiD /2ttT 
is the thermal diffusion length. It has the meaning of the dimensionless conductance of the wire segment of length 
Lt- The other one is the number of transverse channels in the wire, N c h = kpA/An. We consider a multichannel 
metallic wire, for which both parameters are large. From now on the Planck's constant H is set to unity. 

For large G(Lt) we can evaluate the replicated partition function in the saddle point approximation. In this 
approximation the partition function is written as a sum of the contributions arising from all the saddle points: 

( Z ' P ) = e~ s ' p - Jv[SQ,SV]e- ss ^ SQ ' sv ^. (6) 

saddle 



Here S s . p , denotes the NLcrM action evaluated at the saddle point, and 8Q, SV describe fluctuations of the Q-matrix 
and electric potentials V a around a particular saddle point. Finally 8S[6Q, SV] denotes the action change due to these 
fluctuations. In the next section we discuss the saddle points of the action J2J in the G{Lt) = const, N c h — * oo limit, 
which is referred to below as the "Nch —> o° limit" for brevity. 



III. SADDLE POINTS IN THE N ch -> oo LIMIT 

If the number of channels in the wire is sufficiently large, e 2 vA 3> 1, one may neglect the Coulomb action, Sc, in 
Eq. when looking for the saddle points. This corresponds to the charge neutrality limits, which can be seen by 
noting that formally such procedure corresponds to the limit e — * oo, which clearly enforces electroneutrality. The 
saddle point equations in this limit are obtained by minimizing Sq in Eq. (|2a|) with respect to V a and Q, and read 

DX7(QX7Q) -[i + V,Q] = 0, (7a) 
V a -^trQ a T a T (x)=0, (7b) 

where tr is the trace in the S®T space only. Equation (|7a|) is the Usadel equation, and Eq. (|7bj) represents the charge 
neutrality condition. 

By direct substitution one can check that Eqs. ((7|) possess a set of stationary spatially uniform solutions, Q a e \, = 
8 ab 8 es i ToSgn(e + 2irTw a ), V a — 2irTw a , which are characterized by a set of integer winding numbers in each replica, 
w a . All these solutions represent degenerate minima of the action (|2bj) . The sum 4^ a w a = W (the factor 4 here 
arises from the 4x4 matrix structure of Q a e h e , in the S T space) defines the trace of the Q-matrix, TrQ = 2W. 
The Q-matrices corresponding to the minima with different w a , but the same W can be transformed into each other 
via continuous rotations in the replica and Matsubara spaces, Eq. ((4]). Therefore, the NLcrM action contains soliton 
minima in which the Q-matrix and the potentials V a smoothly interpolate between their values in different uniform 
minima 15 . Such solitons are similar to those first found in Ref.— . 



A. Single soliton solution 



In this section we find an analytic solution to the saddle point equations |J7J) that correspond to a single soliton. 
To be specific, we construct a soliton that connects the following degenerate minima: Q = A, with all the winding 
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numbers w a = at x — — oo, and Q ab e , — S ab S eE 'T sgn(e + 2irTw a ) at x = oo, with wi, 2 = Tl, all the other w a being 
zero. This corresponds to a gradual change in the electric potential in replicas 1 and 2, Vi 2, from zero at negative 
spatial infinity to =p27rT at positive infinity. 

For such a soliton the generator Wq parameterizing the saddle point Q-matrix via Eq. Q corresponds to a rotation 
between Matsubara frequencies ttT in replica f , and — irT in replica 2. In this subspace Wq has the following structure: 



Wn = 



A 

A+ 



A = 



9 d e* 



6 c e*x 
-Ode'"' 



(8) 



where 9d, 9 Cl <f> and x are rea l parameters. In this equation the matrix element of Wq in the upper-left-corner cor- 
responds to (VF)^ nT = 0, the one in the upper-right-corner to (W)*y T — A, and so on. All the other matrix 
elements of Wq are zero. 

Substituting the rotation generator (|8|) into Eq. ([4]) we obtain the matrix elements of the Q-matrix that participate 
in the rotation: 



O 12 



Q^ttT.ttT II Q Z -ttT,-ttT 



7T 



cos 6d cos 8 C 


e i{<t>+X) sm d S i n c 


—ie 1 ^ sin 6^ cos 9 C 


—ie lx cos 9d sin 6* c 


e -i(4>+x) S i n d s - m c 


cos 6d cos 8 C 


ie~ %x cos 0d sin 9 C 


ie~ 11 ^ sin 9d cos (9 C 


ie~ l<t> sin 9d cos 8 C 


—ie lx cos 0d sin 9 C 


— cos cos 9 C 




ie~ %x cos 9 d sin 9 C 


~ie 1 ^ sin 9d cos 9 C 


e i(4>-x) s i n d s [ u q c 


— cos 8d cos # c 



/ 



(9) 



All the other matrix elements are those of the A-matrix. 

The action for such a Q-matrix is independent of the angles <fi and % anc l depends only on (V0) 2 and (V%) 2 
with positive coefficients. Therefore the action minimum corresponds to coordinate independent angles <f> and \- It 
can be shown that the soliton solutions with the minimum action correspond to either 9d 7^ 0, 9 C = (diffusonlike 
rotation) or 9 d = 0, 8 C 7^ (Cooperon-like rotation). In these cases substitution of Eq. © into (|7b[) gives Vi^(x) = 
=F7rT[l — cos 6d, c(x)] for the diffusonlike and Cooperon-like rotations respectively. Then Eq. (I7ap yields 



1 



2Lrp 



sin 29 d 



0. 



The solution that corresponds to the sought soliton is 

9 d ,c(x) = 2arctan ( e (:Mo)/LT ) = 6> (a; - ar ), 

giving for the electric potentials 

Vi, 2 (x) = tV°(x - x ) = TttT{1 + tanh[(z - xq)/Lt]}, 



(10) 



(11) 



(12) 



which clearly satisfies Vi^(x — > —00) = and Vi,2(x — > 00) = =F27rT. Here xq denotes the soliton position. 

Substituting the saddle point values of Q and V a , Eqs. (J9j> and (fl2|) . into the action (|2bj) . we obtain the action for 
a single soliton 

So = G{L T ). 

We note that this action does not depend on the soliton position xo, and the angles 4> and X m Eq. (|9|). However, for 
the diffusonlike (9 C = 0) soliton the different values of the angle x correspond to the same Q-matrix, and similarly 
different values of <j> correspond to the same Q-matrix for the Cooperon-like (9d = 0) soliton. Therefore the action for 
the fluctuations about the soliton has only two zero modes. One is associated with a translation of the soliton (change 
in xo). The other corresponds to a rotation of the Q-matrix in the replica and Matsubara space caused by a uniform 
change in either 4> or depending on whether we consider a diffusonlike or a Cooperon-like soliton. The presence of 
these zero modes needs to be borne in mind when integrating over the fluctuations about the soliton configurations. 



B. Fluctuations around a single soliton 

In this section we evaluate the single soliton contribution to the replicated partition function, Eq. ([5]), in the 
N c h — > 00 limit. This requires evaluating the functional integral over the fluctuations of the Q-matrix and the 
potentials V a around the single soliton saddle point. 
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As was explained at the end of Sec. IIII A[ the fluctuation spectrum has two zero modes. We show below that all 
the other fluctuations are massive and integrate over them in the gaussian approximation. The resulting fluctuation 
determinant is convergent and is evaluated below. The integration over the zero modes is reduced to the integration 
over the soliton position and the rotation angle. 

The translational zero mode represents a simultaneous spatial shift of the saddle point solution for the Q-matrix 
and the static (zero Matsubara frequency) component of the potentials V a - In order to simplify the treatment of this 
zero mode we first integrate over the latter. This step involves no approximations since the action (|2b[) is quadratic 
in V a . The resulting action depends only on the nonzero Matsubara components of V a and on the Q-matrix. In this 
representations the zero modes involve only the Q-matrix degrees of freedom, whereas all fluctuations of the nonzero 
Matsubara components of V a are massive. Then the single soliton contribution to the partition function, Eq. ^ , in 
the N c h — > oo limit can be written as exp[— G(Lx)] T p , where T p is the functional integral over the fluctuations about 
the soliton and is given by 

r p = a* [v[W,6V]e- s(2) W sv l (13) 

Here the fluctuations of the nonzero Matsubara components of the electric potential are denoted by 8V a , the matrix 
W parameterizes the deviation of the Q- matrix from the saddle point, and a p is the factor coming from integration 
over the static components of V a . We will see later that in order to obtain the physical observables we will only need 
to evaluate Y p at p = 0. Therefore, the value of a is of no importance. Finally, the quadratic fluctuation action 
[W, SV] is obtained by integrating over the fluctuations of the static component of V a in Eq. (|2b[) . and expanding 
the resulting action to the second order in W. Its form depends on the Q-matrix parametrization. 
In the remainder of this section we show that the fluctuation integral T p can be expressed as 

T p = a?G{L T )T p j ^, (14) 

where x$ is the position of the soliton and T p is a numerical factor independent of the system parameters. In order 
to evaluate the thermodynamic quantities we need only the p = value of this quantity, which is calculated below, 
T = T p=0 «8. 

In the remainder of the present section we derive Eq. fT4")) . The presentation is organized as follows. In Sec. IIII B II 
we give the expression for the fluctuation action. In Sec. IIII B 21 we carry out the integration over the Q-matrix 
fluctuations. Section fill B 31 deals with integration over the electric potential fluctuations. The reader not interested 
in the derivation of Eq. (I14| may wish to proceed directly to Sec. IIV| where we use it to evaluate nonperturbative 
corrections to the thermodynamic quantities. 

1. Fluctuation action 

We parameterize the deviations of the Q-matrix from the saddle point in terms of the matrix W, whose structure 
is described by Eq. §5$, as follows: 

Q = e W /2 e W/2 Ae -W/2 e -Wo/2_ (15) 

Here the matrix Wq parameterizes the saddle point Q-matrix. For the soliton described in Sec. IIII Al it is given by 

Here i = 0, 1 corresponds to diffusonlike and Cooperon-like rotations, #o is defined in Eq. (fTTj) , and we set <fi = \ — 7r/2, 
and Xq = Q for convenience. 

In the following we use dimensionless coordinate £ = x/Lt, dimensionless fermionic Matsubara frequencies, e = 
e/2irT, and dimensionless Matsubara components of the electric potential, V" = 5V^ /IttT, where ui is an integer 
defining the bosonic Matsubara frequency, such that the latter is written as 2-kTlo. In these variables the quadratic 
action in Eq. (fT3"|) can be written as 

S {2) [W,V] = s vv + s ww + Swv- (17) 
Here Syv denotes the part of the action that is quadratic in the potentials V a , 

S VV = EE/ W(0V o - w (r), (18a) 



6 



Sww denotes the part of the action that is quadratic in W 
G{L T 



Sww 



16 



(cos 9 — 1) 



/ d q E E {( £ - ^')tr^(^)t + trW^VCW.f,) 1 } 

^(tr^ )£ (W^ e )t + trW e a ?_„ T (W e % T ?) 



ab e>0 : e'<0 
(1 — COS 6*o ) 



]T sgne(tvW e \ a , (W^ - txW% (W^ ) t ) 



sin 2 6»o 



tr {(r,^.^) 2 + (r^ 21 ^ T ) 2 } - - (tr^W^.^ - W^ T )}Y 



(18b) 



and Swv denotes the part of the action that is linear in W and V, 



iwv 



.G(Lq 



. G{Lt) 



e+2-nTuj 

/^E {- ( cos | - i) ^Kt^-^ +sm|tr[r^ 2 12 T(i 



,-ti-tJ 



-Vr w (-sin— trbl^, 12 1-L , J + ( cos — - 1 ) trW" „ nl A 



cos - 1 ) ti\W 2 l rpfi , J f - sin -^tr^W/ 12 , „,,,. , ,1 



2ttT(^+u;)J 



-27rT(i-w),-7rT 



(18c) 



Here "f" denotes the Hermitian conjugation in the S(&T space, i.e. corresponds to complex conjugation and transpo- 
sition within a 4 x 4 block, without interchanging replica or Matsubara indices. The diffusonlike soliton corresponds 
to r,; = tq, and t, = T\ for the Cooperon-like one. To be specific, in what follows we consider the case of a diffuson-likc 
soliton, i.e. we set t\ = to. In the Cooperon-like case the treatment exactly parallels the one presented below. 
Introducing the notation 



T w = J VWeM-Sww), 



and 



iv = J v[v] 



e -S VV / e -Sw 



V[V]e 



where (. . .) w denotes the gaussian average with respect to the action SwWi we can write Eq. (| 13() as 

T p = a p Tw^v- 

We evaluate quantities Tw and IV in Sections IIII B 21 and IIII B 3l 



(19) 



(20) 



(21) 



2. Integration over W 

We now evaluate the functional integral over the fluctuations of the Q-niatrix, IV in Eq. (|19p . Examination of 
the quadratic action in Eq. (|18b|) , shows that the variables W^, with different replica or Matsubara indices fluctuate 
independently. Moreover, with the exception of W^j- _ tT , for each W~l the actions for the diffusons and Cooperons 
constituting it are identical. The term containing _ 7rT is special because it has the same replica and Matsubara 
indices as the rotation generator Wq parameterizing the saddle point. The fluctuations of the diffuson and Cooperon 
components of _ nT are also independent, but their propagators are different. In particular, we will see that for 
a soliton represented by a diffuson-like rotation only diffuson part of T has zero modes, and vice versa for a 

Cooperon-like rotation. 

In terms of the diffuson and Cooperon variables, see Eq. (JSJ), the action (| 18b[) can be written as 
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TABLE I: Potentials u£, appearing in the operators L% e i, Eq. (|23p . Each entry gives the potential specific to particular replica 
and Matsubara indices in terms of the potentials v\,% and u defined in Eq. (|26[) . The Latin letters (j, k) denote replica indices 
not equal to 1 or 2. 



ee' \ ab 


jk 


y 


2j 


Jl 


J'2 


11 


12 


21 


22 


e > nT,e' < -nT 





«2 — 1 


Vi 


W2 


H2 — 1 





2^2 - 2 


2wi 





e = 7rT,e' < -ttT 





?i 


«l 


^2 


^2 — 1 


Wl + it 


t; 2 + u — 1 







e > 7rT,e' = -ttT 





— 1 


«i 


^2 


It 





^2 + u — 1 


2uj 


«i + u 


e = 7rT, e' = —ttT 





It 


Vi 


V2 


It 


fl + u 


Excluded 


2«j 


«1 + u 



26 

e>0,e'<0 



where the primed sum means that the term with a = 1,6 = 2, £ = 7rT and e' = — 7rT is excluded, and (d, c) s 
(<i, c)*y -ttt- The operators -£^, £d,c are all of the Schrodinger type and have the form, 

G(L T ) 



Ld,c = — j — \ L u =i + Ud,c{£) ) , 



G(L T ) 



(ie-e' + C^CO). (23) 



with the operator defined as 



d 2 , , 

L » = -^+"> ( 24 ) 

with ui and e being the appropriate dimensionless Matsubara frequencies. The potentials u ( i. c for d s , c s are given by 

Ud(0 = ~c^n%)' Mx) = ~^-y (25) 

The potentials U^, depend on the replica and Matsubara indices involved and can be expressed in terms of the 
following potentials, 

« 1 ^(0 = i[l±tanh(0] 1 ^) = - 4cQs 3 h2 ■ (26) 

The expressions for the potentials £/? e h , in terms of Vi^iO and u(£) are summarized in TableQ] 

The operators L^, and L c are positive definite. The operator L^, Eq. (j2"5)) , with the potential u^, defined in 
Eq. (|25j) . has one zero eigenvalue, with all the other ones being positive and separated by a finite gap. The integration 
over the zero modes requires a special consideration. We therefore defer the integration over the variables d s in T\y, 
Eq. (|19p . to the end of this section and begin by integrating over all the other variables first. To this end we introduce 
an auxiliary quantity T' w as 

j V[d s ]e- $ dtLd *^ dds 

Tw = ~r — ; — ; g(-l t) ,. c ,, f —T' w = T d T' w . (27) 

Calculation of T' w reduces to evaluation of gaussian integrals. Since (d, c)"|/ are complex fields, the integration over 
each of them gives a factor of an inverse determinant of the corresponding operator in the quadratic action, Eq. (|22j) . 
and we obtain the following expression for T' w , 

p „ / 

ab 



^w = TTrTT 7 \ II ( det ^ 



ab 

e>0,e'<0 
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The prime indicates that the product does not include the contribution from a = 1,6 = 2, e = irT,e' = —ttT. The 
operators L^, in the expression for T' w can be classified according to whether their replica indices correspond to 
the replicas participating in the soliton rotation. In particular, for a, b > 2, the operators L^, are insensitive to the 
presence of a soliton. Denoting each of these operators as L? ,, we see that the product over the replicas with a, b > 2 



contributes a factor I J] 



le>0,e'<0 



det L 



-(p-2Y 



to the fluctuation determinant. Analogously, for a = 1, 2 and b > 2, 



we have p — 2 identical operators L a ^ e , for each of a = 1 and a = 2, which we denote as Lj £ , and L^, respectively. 
Finally, there are p — 2 equal operators for a > 2 and each of 6 = 1 and 6 = 2, denoted as 1? , and Lg S /. Using these 
observations we rewrite the previous equation as 



r', 



a 



v 



\\~ 



det 



£a>=i ) detL c 



\\ det L^, det Lf e , det L^, det Lg, 



4 ^ u =x j ■-'c \ £ >0,e'<0 
-(P-2) 2 



n detL ^ 



n det L« det det det i£ 



, e>0,e'<0 



. e>0,e'<0 



-(p-2) 



In the above expression the prime means that det L X2 T _^ T is excluded from the product. To compute the thermody- 
namic quantities we will need only the value of T' w at p = 0, for which we use the same notation, 



F' — 
1 w — 



<let(^ii a ,=i)d<-tl 



[] det Ll 2 e , det Lf e , det det L 22 , 



detL5 



c \e>0,e'<0 



x I [] det L^, det L^, det L^, det L^ 2 , 

|£>0,E'<0 



(28) 



Using Eq. flSbj) . the definit ions 
for 



, and the identity In det O = tr In O we can write 



lnlV = 2^ (4wtr £ ln 



(L u +Vi)(L u + v 2 ) , . {L u + 2v 1 ){L UJ + 2v 2 ) , . L, 
— wtrj m ~ r h 4tr^ In — 



-2tr c In 



L U (L U + 1) 
(L^ + v 1 + it)(L^ + t; 2 + u) 



tr^ In 



L u (L u + 2) 



(29) 



where tr^ denotes the trace in the coordinate space, tr^O = J d£0(£, £). The terms in Eq. (|2"9"|) are grouped in such a 
way that each is finite at a given to, i.e. does not diverge with the length of the system. 
In Appendix [A] it is shown that each term in Eq. I|29p can be evaluated using the formula 



tr^ In 



[ l L u + U 1 )[ l L U) + U 2 ) 
L U (L U + h) 



= hxtt' = In 



(30) 



where the potentials U\{£) an d U 2 (£) satisfy J7i(£) = U 2 (—£,), U\{— oo) = 0, U\(oo) = h (for the potentials from 
Eq. (|29p the parameter h takes on the values 0, 1 or 2). The quantities t,t' describe the £ — ► oo asymptotics of the 
two independent solutions of the equation 



t/> = 0. 



(31) 



Namely, if we find the two solutions, tpi l2 , whose asymptotics at £ — > ±oo are given by ipi(£, — ¥ —oo) ~ ex p(V^0 
and 1^2 (£ — > +oo) w exp(— V^T+Ti £), the parameters f and i' are given by the coefficients in front of the growing 
exponentials in the asymptotics of these solutions at the opposite infinities, 



V>i(£) ~ texp(VuJ + hg), £-+oo; ip 2 (0 « *' exp(- Vw£), £ -> -oo. 



(32) 
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The last equality in Eq. (|30l) holds since \Juj + ht — yfuit' , see Appendix [XI for details. The case of a potential vanishing 
at spatial infinities is recovered from Eq. (|30|) by setting h = 0, U± = U% = U, t = t': 

tr e ln X " + [/ =lnt. (33) 

In order to find the parameters t and t' corresponding to the potentials in Eq. ([29]) . we note that for each potential 
Eq. (|3~lj) has the general form 



2 



cosh 2 £ 2 



^(1+tanhO 



V; = 0, (34) 



where the values of the parameters a, [3 depend on the specific potential. For example, the potentials it(£) and i>i(£) 
in Eq. (|26)l correspond to a — 3/4, (3 = and a = 0, /3 = 1 respectively. 

If one introduces the variable z = (l + tanh£)/2, and y(z) = z~^l 2 {l — z)^^ u)+ ^^ 2 il){z), the above equation reduces 
to the hypergeometric equation for y(z): 

z(l - + [c - (a + b + l)z] ^ - aby = 0, (35) 

where the parameters a, &, and c are given by the following expressions: 

c = 1 + v 7 ^, 
1 



a = -(1 + + + P - Vl + 4a), 



b = -(1 + + y/u + (3 + Vl + 4a). (36) 

Using the properties of the hypergeometric functions F(a, b, c, z)2£ and switching back to the original variable 
£ = arctanh(2z — 1) it is easy to show that the two independent solutions -01,2 of Eq. (|34[) satisfying the desired 
asymptotics, 0i(£ — * — oo) — > cxp( v / a;£) and 02 (£ - * °°) - * cxp(— v 7 ^ + are given by 

0i(z) = z^(l - z)^ E F{a,b,c,z), 

fa(z) = z^(l- z)^ E F(a,b,a + b-c + 1,1- z). (37) 
The asymptotic behavior of tpiiO a ^ £ ^ +°° 1S 

M ^^)« T ^ + l- e K ^, (38) 
T(a)T(b) 

where T(x) is the Euler gamma function. Comparing Eq. (f3"8")) with Eq. (f3"2"|) . wc find that the value of the coefficient 
t entering Eq. (|30[) is given by 

T(c)T(a + b - c) , x 

'- r(V) ■ (39) 

with a, b, c defined in Eq. ((351). Using Eqs. §5U$), (J33J) and and the identity T(x + 1) = aT^x) we obtain for T'^, 
Eq. G2D: 



lnT^ = 2]T Ucjhx 



-2 In 



In r 



4y^yZ7TT r 2 (yzj)r 2 (VZ7TT) 
(•v/w + V^TT) 2 r 4 (v^/2 + V^TT/2) 

i6yW^TT r 2 (yzj)r 2 (V^"+T) 

((V^+V^TT) 2 - 

/ 3- V5\ ^ / 3 



w In 



(V« + V^Tf) 2 r 4 (v^/2 + v^T2/2)_ 



1) 2 r 4 (v^/2 + - 1/2). 

In 0.08. 



41n 



(a; -1/4) T 2 (V^-l/2) 



(40) 
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One can check that the sum over u> above converges, since the summand behaves like uj~ 3 ^ 2 for large u>. The last 
equality was obtained by performing the summation numerically. 

We now complete the evaluation of rv by computing the functional integral Yd, defined in Eq. (|27|) . The operator Ld 
defined by Eqs. (|23[) and (|^5|) has one zero eigenvalue. The corresponding eigenfunction is 1/ cosh£. The fluctuations 
of ^td s and Qd s in the numerator of Eq. (|2T[) along this mode correspond to the rotational and translational zero 
modes of the soliton discussed at the end of Sec. IIII Al 

Indeed, in the parametrization (|15p a soliton displacement, Qo(£) — > Qo(C ~ £o)j by a small amount £o — %o/Lt is 
described by the generator W^ that can be obtained from the condition 

S Q « -i^o/2 AW(oe -iWo/2 = ^ 

where Qo is given by Eq. (fT5|) with W = and Wo from Eq. (fT6|) with i = 0. From this equation it follows that W^ 
has the same structure as Wo, Eq. ((SJ), with matrix A replaced by (W^ )^ T _ 7rT defined as 

( W io)nT,-nT = ~ lT0 —d£-t°- 

Comparing this expression with Eq. we see that the soliton translation corresponds to the diffuson fluctuation 
of the form d s (£) = Co = — (*/ cosh£) £q- Along the same lines of reasoning it can be shown that the soliton 

rotation by the angle <f>o, <fi — > ir/2 + <fi in Eq. corresponds to d s (£) = (1/ cosh£) (j> , and represents the other 
zero mode. 

We separate the functional integral over d s in the numerator of Eq. (|2T|) into a product of integrals over the zero 
and massive modes: 



J J dtj J dfo J g[4]e^fg / dtd s '(L^ 1+ u d )d 3 

where d s contains the massive modes only, and J denotes the Jacobian for the change of variables {d s } — > {d s , £q, 
The product of the Jacobian J and the ratio of the functional integrals in this expression can be evaluated using 
the following trick. We introduce a regularized ratio Yd(rj) of the functional integrals over d s in the last equation by 
infinitcsimally shifting the frequency u> from unity, uj — > 1 + 77, where rj is positive. As a result the zero modes acquire 
a finite mass and 1^(77) can be written as 

r d(v) = L G( , T) f J,. f ; • (41) 



On the other hand this ratio of gaussian integrals can be calculated using Eq. (|33[) . (j3T))) . and In the limit of 

?]->0we obtain 



WJ L„=i+„ + « d r(VT+^)r(VT+^+i) t? 1 J 



To arrive at this expression set u = 1 + 77 in Eq. (|36[) and used the fact that the potential Ud corresponds to a = 2 
and (3 = 0. Integrating over £0 an d <^o m Eq. Pdj) and comparing the result with Eq. (|l2")) we conclude that Td can 
be written as 

Substituting this expression into Eq. (|27p and integrating over <po we obtain the following expression for Yw, 

r w =4T' W G(L T ) J d£ , (43) 

with T' w given by Eq. (TU)) . 
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3. Integration over the electric potential fluctuations 



We now turn to the evaluation of the functional integral over the potential fluctuations, Ty in Eq. (|20|) . The action 
for the potential fluctuations is obtained by evaluating the gaussian average, {Swy) w , in Eq. (12T))) with respect to the 
action Sww hi Eq. (|18b[) . The result of this tedious, but straightforward calculation can be expressed in the form 



G(L 7 



a=l,2 



(44) 



Here, in the p — 2 replicas not participating in the soliton rotation, the dimensionless polarization operator IIq (£ — £') 
is given by the usual expression, 



and in the remaining two replicas (a = 1,2) the dimensionless polarization operators acquire a correction (HT^ (£,£') 
due to the presence of a soliton, 



Sn»(t,e) = cos^G-^Ocos^+sin^G-^nsin^-G^e-O- (46) 
In the last equation we introduced the following Green's functions: 

Gi,a(£>0 = (L M +vi,2 + u)-\ (47) 

where the operator L w and the potentials i>i(£), ^2(£)> an d u (£) are defined in Eqs. ([24)) and (|26]) . 

We note that the polarization operator in the presence of the soliton, Eq. (|4"4")l . is diagonal in Matsubara frequencies. 
This is a consequence of the fact that the soliton saddle point is static. We also note that no inter-replica couplings 
between the potential fluctuations are generated. 

As an important consistency check, let us prove that 

/oo 
dC , n^°(e,0 = o, (48) 
-oo 

which must hold due to particle number conservation. The polarization operator IIq automatically satisfies this 
property, as its Fourier transform is proportional to q. To prove that SIL" satisfies the same condition, we note that 
(—^2 + V\ + u) cos ™ = 0, and (—^2 + v 2 + u) sin ^ = 0. Thus we can write 



f" dgSW (£, O = r d? cos ^ Gt % O cos ^ + j°° # sin ^ G£ (£, £') sin 

J— 00 J— 00 * * </ — 00 * 



- O = n cos 2 I + r\ sin 2 I - ± = 0, (49) 
|w| 2 |w| 2 |w| 

as expected. 

Performing the gaussian integral over V in Eq. (|20j) and taking the number of replicas p to zero we obtain 



IV= n det2n » =e -2E a)>0 tr«in[i-5n-(ns')- 1 ] f50) 

Due to the complicated form of JII" evaluation of this quantity explicitly is a daunting task. In particular, the method 
of the previous section does not apply here because the polarization operators are not represented by Schrodinger-type 
operators. However, we note that the dimensionless polarization operators and — SIl^ are independent of the 
system parameters. Provided the sum over u in Eq. (|50j) converges, it is clear that Ty is a parameter-independent 
numerical factor. Below we prove that the sum over u> in the exponent of Eq. (|50[) does converge, and evaluate Ty 
numerically. 

To this end we obtain the large-w asymptotics of the summand in Eq. (|50[) . This can be done by expanding the 
logarithm in Eq. fl50|) to the first order in iflP^nft) -1 , 
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t^lnfl-OT'W)- 1 ] M-trtSIFQE$)- 1 = -J d^dC'SU^te)^)^. (51) 
Using the Fourier transform of IIq from Eq. (|45[) . the above trace can be written as 

tr^IFfffi')- 1 = f — f dtdge-^SnPte^e**' . (52) 
J 27r q 1 J 

We note that each of the two terms in the expression for (511", Eq. (|46|) . can be written as ipo(QG(£,(i')il)o{£'), 
where C?(£, £') is the resolvent of of the operator (w — ^ + U), and V'o is the zero mode of (-^ + U). The 
phase factors in the last integral in Eq. (|52|) can be interpreted as a gauge transformation of the Green's function 
G(£,g) = e- 4 ««G(£,£')e l9? ' =(u + (i4 + g) 2 + J/)" 1 . Therefore, the integral can be written as 



df^'e^OT'' = J d£d? cob 



9o(0 



w + <7 



2 ' 



1 



, e . #o(£) 
"que; sm — 7, — 



+ q 2 + 2 q C-±)-^+v 1 +u / 

1 



cos 



<? (r) 



. wo 

sm — - — 



(53) 



dc 5 



+ ^2 + 11 + 29(11) /(w- 



Then we expand each of the first two kernels in the r.h.s. of Eq. (|53[) in powers of 

to the order that gives first nonvanishing contribution to the entire integral. The zeroth order term vanishes in the same 
way as it happened in Eq. (|4U)) , and so does the first order one being proportional to cos ^ ^ cos ^ + sin ^ ^ sin ^ = 

4^(1/2) = 0. Therefore, expansion to the second order gives the first nonzero contribution, and we arrive at 



trtSU (K + U ) « y ^^(- 4 )(^Fi ^\ C ° S 2-^ COS 2- +Sm 7^ Sln T/ 



' (54) 



2w 3 / 2 



Eq. (154|) proves convergence of the sum over frequencies in Eq. (|50[) . Therefore we do not have to introduce any 
additional regulators. 

We can calculate IV numerically by expanding the logarithm in Eq. (J5UJ) v in <5IU(r±g) _1 and calculating the 
corresponding traces. The explicit calculation shows that expansion to the third order yields a precision better than 
a percent, which is sufficient for our purposes. Proceeding this way we obtain IV ~ 24. Combining IV with IV, 
expressed via T' w and r<j calculated in Eqs. (|40|) and (|43|) . we obtain the final expression for T, determining the 
fluctuation integral r p= 0i Eq- (|14p . needed to calculate the thermodynamics quantities: 

T = 4? W T V « 8. (55) 



IV. NONPERTURBATIVE CORRECTIONS TO THE THERMODYNAMIC QUANTITIES 

In the previous section we found the soliton saddle points and showed that the functional integral over the fluctua- 
tions about a single soliton configuration can be expressed in terms of the integral over the soliton position, Eq. (|14[) . 
In the present section we use these results to obtain nonperturbative corrections to the thermodynamic quantities at 
relatively high temperatures, G(Lt) 3> 1. We begin by considering the N c h — > oo limit in Sec. II V Al and turn to the 
case of large, but finite N c h in Sec. IIVB1 
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Infinite channel number 



In the N c h — * oo limit the Coulomb action (|2c|) vanishes. In this case the NLcrM action has infinitely many 
degenerate saddle points with spatially uniform potentials characterized by the winding numbers w a , V a (x) — 2irTw a , 
with the usual saddle point, Q = A, corresponding to all w a — 0. The single soliton solutions, studied in Sec. IIII Al 
represent exact inhomogeneous saddle points with a finite action G(Lt) and correspond to a kink- like change of of the 
electric potentials V a (x) by ± 2irT, Eq. (fT2"|) , in two of the replicas involved in the soliton rotation. The spatial size of 
the kink is given by the thermal diffusion length Lt- In the dilute soliton gas limit, which corresponds to G(Lt) 1, 
multi-soliton saddle points can be viewed as sets of such kinks separated by distances much larger than Lt- In this 
case the action of a multi-soliton saddle point is given by the sum of single soliton actions. Similarly, the functional 
integral over the massive modes factorizes into a product of fluctuational determinants for each soliton. Thus in the 
dilute regime the soliton gas is noninteracting. Noting that the sum over the saddle points in Eq. (J6j> factorizes into 
a product of a sum over the uniform saddle points and the sum over the soliton configurations, we can easily find the 
multi-soliton contributions to the replicated partition function in the dilute soliton gas regime, 

(zn = z> ± [^-WMe-°w] B ^ , ^ = zi ^ ^ _ 1KTpG(iT)e _ G(ir)L/iT ] (56) 

n=0 ' i=l ^ 

Here Zfi denotes the contribution of the homogeneous saddle points to the replicated partition function, L / Lt is the 
dimensionless wire length, the factor of p(p — 1) arises from the number of ways the two replicas participating in the 

(i) 

soliton rotation can be chosen from the p replicas available, £q denotes the position of the i-th soliton, and the factor 
of two arises from taking the Cooperon-like and diffusonlike solitons into account. Substituting this result into Eq. ([1]) 
we obtain the leading nonperturbative correction to the average thermodynamic potential in the N c h ~ > °o limit: 

«)«, = 2TG(L T )e- G(LT) ^T, (57) 

Lt 

where T is defined in Eq. (|55|) . 

The correction to the heat capacity can be obtained as <5Coo = — T d . Taking into account that the largest 
contribution comes from differentiating G{Lt) in the exponential, we obtain the ratio of SCoo to the heat capacity of 
noninteracting electrons, Cq = ^-vATL, 

SCry, 



Co 



-24TG 2 (L T ) e - G(LT) . (58) 



The analysis above was restricted to the charge neutrality limit, N c h — ► oo. In the next section we consider the 
case of a large, but finite number of channels in the wire. In this case the Coulomb action (l2"cl) may not be neglected. 
Its presence significantly modifies the behavior of the soliton gas. 

B. Finite channel number 

For N c h ^> 1 the influence of the Coulomb action (|2"c)) on the soliton shape and on the massive fluctuations 
about the multi-soliton configurations is small and may be neglected. Therefore, each soliton configuration is still 
fully characterized by the soliton positions and the indices of the replicas participating in the soliton rotation. The 
Coulomb action for each configuration is given by the term Sc, Eq. (|2c[) . evaluated for the specific potential profile 
V a (x) corresponding to such a configuration. 

For a single soliton situated at Xq the potential profile in the two replicas participating in the rotation is represented 
by a kink, Vq(x — xq), Eq. (fl2|) . in one of the replicas and an antikink, —Vq(x — Xq), in the other one. Thus each 

(i) . . . . (i) 

soliton is characterized by its position x Q and the index of the replica containing the kink, a + , and the antikink, 
. The potential profile for each soliton configuration is given by 



W = E K,a« - * a , a u] V °( X - ^ ( 59 ) 
i 

Using this representation, the replicated partition function, Eq. ([S]), can be written as 

w ,^ Wi- a /«S°e-p[-*({<M})], (60, 

n " a i *f fA\ 



n=0 i=l" U) 
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where Sc ^£,q \ a ±^}) denotes the Coulomb action (|2"c)) evaluated for a given soliton configuration a± H. Since 

the Coulomb action diverges for any uniform saddle point with w a 7^ 0, such saddle point are forbidden, and the 
factor Zq , arising from the uniform saddle points, contains the contribution only from the usual saddle point, Q = A, 
{w a = 0}. 

Equation (j6"0|) is valid in the dilute soliton gas regime, G(Lt) ^S> 1, in which the typical inter-soliton distance exceeds 
the thermal diffusion length Lt- In the following we assume that at these length scales the Coulomb interaction is 

d 2 

screened due to the presence of a nearby gate, so that its Fourier transform is given by K(q) w 1/e 2 ln(^f), where d g 
is of the order of the distance to the gate. Then defining the kink density p a (0 in replica a, 

Pa(0-T, S ^-^)\S a ^-S n ^}, (61) 



+ 



we can express the Coulomb action in Eq. (j6"0)) in the dilute gas limit as 

Sc (K' "*}) - -i^flr ? / «>"«w?>i« - a m 

Equations (|6"0|). ([BT]). and (f6"2")l express the replicated partition function of the disordered wire as a partition function 
of a one dimensional replicated neutral gas of kinks and antikinks interacting via a linear potential. Importantly, the 
positive and negative charges in this gas occur only in pairs, such that the appearance of a positive charge in one 
replica is accompanied by the appearance of a negative charge in a different replica at the same spatial position. This 
problem can be mapped onto a one-dimensional replicated sine-Gordon model^. Below we will not use this mapping, 
but work in the replicated kink gas representation. 

Only the soliton configurations that correspond to a neutral kink gas in each replica give a nonvanishing contribution 
to the partition functions because all non-neutral configurations possess an infinite Coulomb action. The density of 
the kink gas is controlled by the fugasity, YpG(LT)e~ G ( LT K Depending on its value the kink gas can be in two 
different regimes. At high temperatures, for G{Lt) ^ hiN c h, the gas is dimerized. In other words the kinks within 
each replica form a dilute gas of bound pairs of a kink and an antikink. At lower temperatures, ln(iV c ^) > G(Lt) ^ 1, 
the kink pairs overlap and form an ionized plasma. The dilute soliton gas approximation used to derive Eq. (|60[) is 
valid in both of these cases. We restrict our analysis below to the high temperature regime, G{Lt) > \nN c h- 

For G(Lt) 3> 1, Eq. (|60p may be viewed as an expansion of the replicated partition function in the powers of the 
fugasity, T p G{L T )e~ G ^ T \ In the presence of the Coulomb action, the single soliton contribution to the partition 
function vanishes, since the corresponding Coulomb action is infinite. Therefore the leading term in this expansion is 
given by the contribution of two solitons which corresponds to two kink-antikink pairs in different replicas. We shall 
refer to this object as a dimer. 

To evaluate the contribution of a single dimer into the replicated partition function we express the Coulomb action 
(f6"2")l in terms of the kink-antikink separation within the dimer, £ re i, and substitute the result into Eq. (|6T)|) . Denoting 
the dimer center of mass coordinate by £ cm , summing over all possible pairs of replicas that can accommodate the 
dimer, and recalling that each soliton can be either Cooperon-like or diffusonlike we obtain 

(ZP) = zrll + Mp-l)u p r 2 p G 2 (L T )e- 2G ^- j^ d^ cm J ^ rei e-I^IW£* j , (63) 



where we introduced the notation 



J N 



s 2 ln^ N ch 
ttvf G(Lt 



L T , (64) 



that has the meaning of the typical kink-antikink separation within each dimer. Since Ln ~ Lt / yTrli, where r e ; 
is the elastic mean free time, this length scale is much larger than Lt within the validity domain of the NLerM 
description. Therefore, the dilute soliton gas approximation is justified. Performing the integrals over £ cm and £ re ; in 
Eq. JM]) we get 



(Z p ) = Zl (\ + MP - l)^T^G 2 (I T )e- 2G ^)) 



(65) 



This expression shows that the single dimer contribution to the partition function diverges as the length of the wire 
L goes to infinity. From the second term we infer that the spatial density of dimers is ~ jjfG 2 (LT)e~ 2G ( LT \ In the 
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regime G(Lt) ^ InA^ this density is smaller than 1/ijv, and multisoliton configurations appear as a dilute gas of 
dimcrs. 

Since the dimers in the dilute limit do not interact, the integration over all dimer configurations results in expo- 
nentiation of the correction arising from a single dimer, second term in Eq. (|65[) . 



(Z p )=Z P o 



,71=0 



1 



Ap{p-l)a^l^fG\L T )e-^ 



ryp 



4 P (p-i) a p rl 



-G 2 (L T )e- 2GI - L T) 



(66) 



Using Eq. (fT]) and the definition 
thermodynamic potential: 



we get the expression for the leading nonperturbative correction for the 
32, 



L 



1 T 2 - In ^ N ch G(L T )e~ 2G ^> j—T, 



(67) 



where T is defined in Eq. (|55|) . Using this expression we obtain the ratio of the nonperturbative correction to the 
heat capacity to that of noninteracting electrons, 



5 £ = -^^N ch G\L T )e-^ 
Oq ir vp a 



(68) 



Equations (|6T[) and ([68]) are the main results of this paper. These results are drastically different from the 
expressions (1571) and l|58p obtained by taking the formal N c h — > oo limit. We note that the corrections for the 
thermodynamic potential for infinite and finite N c h, Eqs. ([57)) and (I67p . become of the same order at G{Lt) ~ hiN c h, 
when the dimer gas crosses over into the ionized regime. 



V. SUMMARY 



We studied nonperturbative interaction corrections to the thermodynamic quantities of a multichannel disordered 
wire. Within the replica NLcM formalism these corrections arise from soliton saddle points of the NLctM action. 
In the limit of infinite number of channels, N c h, in the wire we obtained the exact single soliton solution of the 
saddle point equations and evaluated the function integral over the fluctuation about the soliton configuration. We 
showed that for G{Lt) S> 1 and N c h 3> 1 nonperturbative corrections to the thermodynamic quantities of the 
system are described by a partition function for a dilute gas of solitons. The latter is equivalent to the partition 
function for a replicated classical one dimensional Coulomb gas. As the temperature is lowered, this gas undergoes 
a crossover from the dimerized regime of neutral soliton pairs at G(Lp) > lnJV c /, to the regime of ionized plasma 
for G(Lt) < In N c h- The crossover G(Lt) ~ lnN c h 1 occurs at temperatures that are parametrically larger than 
those corresponding to the transition from weak to strong localization, G{Lt) ~ 1- This enables one to study this 
crossover separately from the perturbative effects. We specialized to the high temperature regime, G(Lt) ^ ha.N c h 
and obtained the leading nonperturbative correction to the specific heat (relative to that of noninteracting electrons), 
SC/Cq ~ N c t l G 2 (LT)e~ 2G ( LT \ Eq. ((68)) . We would like to emphasize that this correction is drastically different from 
the result obtained by taking the formal limit N ch -> oo, Eq. ((551), SC^/Cn ~ G 2 {L T )e~ G{ - LT) . 

Although our treatment was specialized to the symplectic ensemble, we believe that the mapping of the nonper- 
turbative corrections to the soliton gas and to the replicated Coulomb gas, described by Eqs. (fBT))) holds for all three 
ensembles. Indeed, the existence of soliton minima is generic for all three ensembles^. The mapping to the replicated 
classical Coulomb gas relies only on the fact that the functional integral over the fluctuations about a single soliton 
configuration can be reduced to the integral over the soliton position, Eq. (fH)) . This, in turn, is a consequence of the 
fact that the integral over the massive modes converges, which we expect to be true for all ensembles. 

The generalization of our formalism to the treatment of nonperturbative corrections to the transport characteristic 
is left for future work. 
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APPENDIX A: DERIVATION OF EQ. (f30|) 

In Sec. lIIIB2l we encountered expressions containing determinants of Schrodinger-type operators of the form 



In Do — tr^ In ■ 



In Dh = tr j In 



c ' 2 + Ut) (u 



(Ala) 



(Alb) 



where E/(£) is a potential that vanishes at £ — > ±00, and Ui t 2{0 are step-like potentials, satisfying U\(£) = Uz{— 0, 
[^(-oo) = 0, L/^oo) =ft. 

As explained in the text above Eq. (I33|) . the trace in (lAlaj) can be obtained as a particular case of that in (lAlbl) . 
Therefore, we concentrate our attention on the latter. We first rewrite Eq. (|Alb|) as 



In A, 



trf In 



trf In 



w-^ + ^i-e) 



tr^ In 



)(' 



3? 



d 2 

3^ 



. (A2) 



where O(0 is the step function. The third term does not depend on the potentials, and can be calculated explicitly, 
which is done at the end of this Appendix. The first two terms are equal since Ui{£) — L^C - 0' We denote each of 
them as lnD^i, and proceed to calculate this quantity. 
To compute 



lnD/ji = tr^ In 



(A3) 



we represent the potential asasum — hQ(^)-\-v(^), where v(£) vanishes at spatial infinities, and express the 

(A4) 



variational derivative of In Dhi with respect to u(0 in terms of the Green's function G(£, £') = — ^ + U\ 
The Green's function G(£, £') can be found by solving the differential equation 



(A5) 



with the boundary conditions that G(£, £') vanishes at spatial infinities — > ±00. It can be expressed^ 3 - in terms 
of the two independent solutions of the homogeneous equation 



mo = 0, 



such that (pi (£ — > —00) — * and (f>2(£ — * +00) — * 0. In particular, at coinciding points we have 
where W[cj)i(^), 02(0] ^ s the Wronskian of 0i(O and 02(C)? 



(A6) 



(A7) 



(A8) 



The Wronskian of the two independent solutions of Eq. (|A6[) does not depend on coordinate £, and therefore may be 
expressed in terms of the £ — > ±00 asymptotics of 0,(0- By appropriately normalizing the solutions, we can express 
the latter as 



0i (0 = 



—00 
00 



, <M0 



i'e- fc? , £ -> -00 



(A9) 
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where k = y/co, c = \Juj + h, and t, t' depend on the specific form of the operator in Eq. (IA6|) . Evaluating the 
Wronskian at £ — > ±00 using the asymptotics (|A9p . we obtain 

Wi(0,^(0] = 2c*=2fc*'. (A10) 

Next we prove that 

S\nD hl Shit 



Sv(0 6v(tY 

To this end we introduce an auxiliary construction 



(All) 



wiMOMO] - - MO^P- (A12) 

Here <j>i and <f>2 are solutions of Eq. (|A6|) with the same w, but for two different potentials and both of which 
vanish at £ — > ±00. The tilde denotes quantities corresponding to v. We assume that <f>\ and 4>2 have the asymptotic 
form (|A9[) , with ^2 characterized by P. 

In contrast to the Wronskian W[^i(£), <fe(OL built out of the solutions of the same equation, the quantity 
W[4>i(£)t 4>2{0] depends on the coordinate and satisfies the differential equation 

d £ = dS 2 - MO dS 2 = K0 - nOJ <p\<P2, (Ai3) 



that follows directly from Eqs. (|A6|) for <j>\ and 02- 

Integrating Eq. (|A13|) with respect to £ from —00 to 00 and using the asymptotic form of 4>\ and <f>2, Eq. (|A9|) . we 
obtain 

/>oo 

W(oo) - W(-oo) = / d£(v - v)(j)i4>2 = 2<rf - 2kt' . (A14) 



Taking a variational derivative of this equation with respect to v(£) at = we obtain 

fc(0fc(O=2*L (A15) 

Plugging Eqs. (|A10|) and (|A15|) into Eq. (|A7f for the Green's function, and using Eq. (|A"IJ) . we obtain Eq. (lAlip . 
Integrating Eq. (| Al 1 [) with respect to v from =0 to its final value we obtain 

hxD hl = hi^, (A16) 
to 

where to is the coefficient in front of in asymptotic form (|A9[) of </>i for = 0. The latter can be easily found 
from the continuity of the logarithmic derivative dln0i(£)/d£ at £ = 0, and is given by to = (1 + k/<;)/2. 

Finally, the third term in Eq. (|A2[) can be calculated in the following manner. We denote this term by T^uj), and 
introduce the Green's functions go, gn, and g that vanish at £, £' — * ±00 and satisfy the equations 



^ 9o = *(S - C), (u-ir + h)g h = 8{S- £'), 



' /2 + he] g+ = 6(Z - £'), ( u -* +h (l-e))g- = 5{£ - £')• (A17) 



Taking the derivative of 13(0;) with respect to u> we obtain 



af = y_ de[ff + «,o+5-(e,o-sbK,o-fffc(e.o]- (ais) 

All Green's functions entering this equation are easily calculated using the method of Wronskian, as was done above 
for a general potential. Specifically, we obtain the following expressions for the Green's function at coinciding points: 



(A19) 
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Keeping in mind that T 3 (uj — > oo) — > 0, we can express it as T 3 (uj) — — duo' dT Q^ - ■ Substituting the expres- 
sions (|A17|) for the Green's functions into Eq. (|A18j) . we obtain 

^3=ln^, (A20) 

with to defined below Eq. (|A16j) . 

Substituting Eqs. (|A20j) and (|Al6jl into Eq. (|A2|) we obtain the final expression for trace in Eq. (|Alb|) . 

- -tm + Ui) ( lu - -£r + U 2 ) fTTTh 
tr* In J \ . * . ' = 2 hiD hl + T 3 = In J ^-^t 2 = In tt' , (A21) 



where in the last expression we used yu + ht = \Zuit' to write a more symmetric expression, in which t and t' are 
defined in |(A"9]) . This proves Eq. (|30|) . 
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